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We have analysed pion photoproduction imposing constraints from fixed t dispersion relations 
and unitarity. Coupled integral equations for the S and P wave multipoles were derived from the 
dispersion relations and solved by the method of Omnes and Muskhelishvili. The free parameters 
were determined by a fit to the most recent data for tv^ and tt*^ production on the proton as well 
as 7r~ production on the neutron, in the energy range 160 MeV < < 420 MeV. The lack of high 
precision data on the neutron and of polarization observables leads to some limitations of our results. 
. . . Especially the multipole Mi_ connected with the Roper resonance Pii(1440) cannot be determined 

' to the required precision. Our predictions for the threshold amplitudes are in good agreement with 

both the data and chiral perturbation theory. In the region of the A(1232) we have determined 
the ratio of electric quadrupole and magnetic dipole excitation. The position of the resonance pole 
is obtained in excellent agreement with pion-nucleon scattering, and the complex residues of the 
^ |i multipoles are determined with the speed-plot technique. 
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I. INTRODUCTION 



Pion photoproduction was extensively studied both theoretically and experimentally in the sixties and seventies. 
The basic theoretical approach was given by dispersion theory, especially by fixed t dispersion relations. Since 

■ perturbative methods of quantum field theory fail in the case of the strong interactions in the resonance region, 
' dispersion relations seemed to be the only possible tool to treat hadronic processes on a quantitative level. Given the 
: technical possibilities of that time, the statistics of the data was moderate and systematical errors were large. With 

On ' the development of quantum chromodynamics (QCD) as a fundamental field theory in terms of quarks and gluons, the 
(-H , theoretical understanding of strong interactions has improved at the higher energies. However our understanding of 

■ hadron physics in the nonpcrturbative region at energies of the order of 1 GeV remains still unsatisfactory. Although 
quark models of different kinds give some description of the baryon and meson spectra, their predictive power for 
reactions is still quite poor. The only approach to infer all the symmetries of QCD into the physics at these energies 
is given by chiral perturbation theory (ChPT). However, the results of ChPT can be applied only in the threshold 
region. Recently, it has been investigated how to include the most prominent resonance of the nucleon, the A(1232), 
within the framework of heavy baryon ChPT. However, the first results achieved by this technique are not yet 
satisfactory. On the other side modern electron accelerators have provided a host of high precision data in the 90's, 

5J] and will continue to do so. Beams of high current and high duty factor together with considerably improved particle 
\ detection techniques have reduced the statistical errors to the order of a few percent, and promise to keep control of 
the systematical errors at the same level. To interpret these data with respect to the most interesting features, i.e. 
the threshold behaviour and the electromagnetic excitation of resonances, a partial wave analysis is mandatory. To 
ensure the consistency and uniqueness of such an analysis, constraints from unitarity and dispersion relations have to 
be imposed. Such concepts have proven to be quite successful in pion-nucleon scattering In comparison with that 
field, the situation in pion photoproduction is considerably more complex. The spin and isospin structure leads to 
twelve independent amplitudes, while in pion-nucleon scattering there are only four such amplitudes. As a consequence 
a complete experiment requires the use of many polarization observables. Such a complete experiment has not 
yet been performed. However, the new experiments provide an ever increasing amount of precise and new data. At 
present, the experimental thrust is mainly on measurements near threshold and around the A(1232) resonance. In 
the coming years, a series of experiments at Jefferson Lab will cover the whole resonance region. Restricting our 
theoretical investigations to the threshold region and the low-lying resonances, we are lead to choose the method of 
Omnes and Muskhelishvili to analyse the existing data, because it introduces a natural parametrization and fulfills 
the constraints of unitarity at the same trace. 

In the future, the required shift to include also higher energies will make it necessary to take account of the full 
content of dispersion methods and to apply more involved techniques. At the moment, however, we obtain an excellent 
representation of the analysed data by use of the method of Ref . M . Therefore we deem it worthwhile extending this 
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method to the second resonance region as a next step. 

In Sec. 2 of this contribution we briefly recall the kinematics, amplitudes and (polarization) observables of pion 
photoproduction. The essential ingredients of dispersion relations at fixed t are outlined in Sec. 3. The method of 
Omnes and Muskhelishvili to solve these coupled integral equations is given in Sec. 4. In particular, we also discuss 
our choice of the free parameters of that method, which arc then fitted to our data basis (Sec. 5). This fit defines 
the result for the multipole amplitudes, which are presented in Sec. 6. In particular we give our prediction for the 
threshold production of charged and neutral pions, the ratio of electric quadrupole and magnetic dipole radiation in 
the region of the A(1232), and the inclusive cross sections. We conclude with a short summary and outlook in Sec. 7. 

II. FORMALISM 

A. Kinematics and conventions 

The four- momenta occuring in the reaction are denoted hy q = {oJ^q) for the photon, k — (uo-^^k) for the pion. 
Pi = (P?, Pi) and Pf = (Pp P/) for the nucleon in the initial and the final state, respectively. With this notation the 
conventional Mandelstam variables are 

s={P, + qf = {Pf + kf (1) 
t=[k-qf = {Pf - P0^ (2) 
U^iPf- qf = (P. - k)\ (3) 

where W is the cm energy. Since all calculations are performed in the cm frame, we further introduce the cm energy 
of the nucleons, 

Eu = Pfr- (4) 
In the following we will suppress the index cm for the momenta, 

g = <f=™--Pr, (5) 
fc = fc"'" = -P;™, (6) 

and introduce the quantity x = cos instead of the scattering angle 9 — 0'^™ . In the discussion of our results for the 
amplitudes and observables, the photon energy is usually given in the lab frame, Ej = EI^^ . 

B. Amplitudes 

The invariant amplitudes Ai[s,t^u) are defined according to The relativistic invariants are|^ 

Mi=i-f^Ell (7) 

M2 = 2i^5{P-qk-e-P-ek-q), (8) 

M3 = ij5{^ k-q-^k-e), (9) 

M4 = 2i75(^P-<Z-^P-e-mAr^^), (10) 

where e is the polarization vector of the photon and P — ^{Pi + Pf). The scattering operator T, to be evaluated 
between the Dirac spinors of the nucleon, is expressed in terms of the invariants by 

4 

T = Y,Ms,t,u)M,. (11) 
1=1 

The photoproduction amplitudes Ai can be decomposed in isospace according to 



^For the metric and the 7 matrices we adopt the conventions of Ref. 
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A, - iA|"^[r„,ro] + A^+^5oO + A|°^a, (12) 

with the Pauh matrices Tq, acting on the isospinor of the nucleon. The behaviour of the invariant ampUtudes under 
crossing is 

Ai{s,t,u) = e'£,kAi{u,t,s), (13) 

/ = +,0,-; ^1 =^2 = -6 =?4 = 1; e+=e" = -£- = l. 
The physical amphtudes may be expressed by the isospin combinations 



A(7p - 


* n7r+) 


A{-fp - 


-.p7r°) 


A['-^n — 




A{'^n - 





(14) 
(15) 
(16) 
(17) 



The amphtudes A'^^^ and A^^) refer to final states of definite isospin (^ or |). They are combinations of the quantities 
AM and A(-) of Eq. (|l|), 

A(^) = A(+)+2A(-), = (18) 

For the isospin i channel it is also useful to define the combinations 

+ „A(^)-A(")-iA(i), (19) 

3 3 

which enter into the amplitudes of the reactions on the proton or the neutron, respectively. 

Evaluating Eq. (|ll|) between the nucleon Dirac spinors in the cm frame, expressing the result by Pauli spinors and 
multiplying by a factor mAr/47riy, we obtain the form 

T — i<7 ■ eJ-i + d ■ k a ■ [q X e)J-2 + id ■ qk ■ sT-^ + id ■ kk ■ eT^, (20) 

where q and k are unit vectors in the direction of the respective momenta. The cm differential cross section can be 
expressed through the CGLN amplitudes JF,; defined by Eq. (|2C|), 

where Xf i ^re nucleon spinors. The connection between the CGLN amplitudes and the invariant amplitudes is given 

by 

2W Ti 



Fi = 47r 



{W-ruN) [{Ef +mN){E,+mNp 



2[W — niN) 



{W -mN) \E^ + mN J \k\ 

f _ 'l7^2 

= -A, + {W + m^)A,~ - A A,~A,), (22) 

2{W + mN) 

i'a = 47r 



Fa = 4tt 



{W-tun) [{Ef +mN){E,+mN)]^ \ k\ 
{W-mM)A2 + {As-Ai), 

2W fEf+niNY -^4 



(PF-TOjv) \E, + mN J I /fc |2 
{W + mjv)A2 + (^3 - Ai). 
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For practical calculations it is useful to introduce a matrix form for the different amplitudes, 



T = llF, F^CA, A^C^'F, (23) 



where A, T and F are vectors with the four elements defined in Eqs. (|n|), ( gO| ) and ( p2| ) respectively. The 4x4 
matrices C and H contain kinematical factors only and are given in appendix ^ 

The CGLN amplitudes are expanded into the multipoles Ei± and M/± which conserve parity and total angular 
momentum J, 



l>0 

+ [{1 + l)Mi^{W) + Ei^{W)]Pl_,{x)}, 

T2{W, 0) = ^[(; + l)Af,+ (T4^) + lMi^{W)]P!{x), 
i>i 

J^s{W,e) = J2i^Ei+{W) - Mi+{W))Pl[,ix) (24) 
+ {EUW) + MUW))Pl'_,i^)l 

Ti{W,e) = J2iMi+iW) - Ei+{W) ~ Mi^iW) - Ei^{W))P;'ix). 



l>2 



E and M denote the electric or magnetic character of the incoming photon, the indices Z± describe the coupling of 
the pion angular momentum I and the nucleon spin to the total angular momentum J = I zt ^. 
The 4 partial waves with angular momentum I may be written in terms of a vector, 

MiiW)^ = {Ei+iW),Ei^iW),Mi+{W),Mi^{W)), (25) 
where T stands for transposed. The matrix form of the multipole expansion Eq. ( [2^ ) and its inverse read 

oo 

J^(M^,a;) =^G,(x)A^i(W^), (26) 
1 

MiiW) = j dxT>i{x)J^{W,x), (27) 
-1 

with the 4x4 matrices G; and D; given in appendix 

The multipole amplitudes AAi± are complex functions of the cm energy W . Below inelastic threshold, the Fermi- 
Watson theorem Q allows one to express the complex phases of the multipoles by the corresponding pion-nucleon 
scattering phase shifts. 



Ml±{W) =\ M\±{W) I eK*i±(w)+"-), (28) 



with I = 



C. Observables 



If polarization degrees of freedom are taken into account, the differential cross section of pion photoproduction 
takes the form Q 

da \ k\ 

— = y^^{{Rt + Pu^jRt"') + ^tHRtt + P„,,,i?TT ) cos2^ (29) 
-(Fi^^^i??^;/ + Pt^^i?^;/') sin2^] + nc(Pi,,^i?TT' + -Pt^./^TT')}- 
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The various response functions R may be expressed through the CGLN amphtudes. Pt are the components of the 
nucleon polarization vector, Tlx is the degree of hnear polarization of the photon, ip the angle between the reaction 
plane and the photon polarization, and He the degree of circular photon polarization. The observables containing 
the information of the various polarization measurements are the appropriate asymmetries. In the following we will 
restrict the discussion to single polarization observables only. The beam asymmetry E which is measured with linearly 
polarized photons, perpendicular and parallel to the production plane is 

S = ^^i^ = (30) 
da-'- + da" Rt 

The target asymmetry T and the recoil polarization P, 

rfa^'H=i+da^".=-i Rt' 

_ da 1 -da f _ Rj, 
da "f~ + da^"f~ ^ Rt 

are determined with respect to the polarization of the incoming (outgoing) nucleon perpendicular to the production 
plane. In appendix ^ the relevant response functions are given as functions of the CGLN amplitudes. 

III. DISPERSION RELATIONS AT FIXED t 

The starting point of our analysis are dispersion relations at fixed t for the invariant amplitudes ||] , 



oo 

ReAiis^t) = 4P°^'(s,t) + / ds' + lmAi{s\t), 

TT J \S' — S S — U J 



(33) 



"thr 



with / = 0, ± and s^j^p = (wat + ttZtt)^. The pole term contributions y^^'P^^^ can be obtained by evaluating the Born 
approximation in pseudoscalar coupling, 

^/,pole,^ ,,_e9f 1 , 



A--- is, t) = ( -— ^ + — ) . — ^ ( ^ + -— ^ ) , (34) 



A 





2 Vs- 


- m- 
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t-ml\ 


i^p°l^.,t) = 


eg 


^1 








eg 




2m N 





2 2 



eg ( e — I e 



1 



ij^ U ^7V 



with K^^"'^^ — Kp — Hn and k^*'-' = Kp + k„, where Kp and k„ are the anomalous magnetic moments of the proton and 
the neutron, respectively. For the pion nucleon coupling constant g we take the value g^/47r = 14.28 |jl[. 

The next step is to apply the multipole projection to the dispersion relations (js^). With the help of Eqs. (|2^) we 
obtain 

oo 

ReF'{W,x) =F''P'^^^{W,x) + ^V J dW'E' {W,W' , x')ImF^ {W , x'), (35) 

^thr 

where x' = cos 9' corresponds to W at fixed t. With 

t-ml = -2luuj^ + 2 I fc 1 1 g I a; = ~2uj'lu'^ + 2 | fc' | | g" | a;' (36) 

we have 
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+ e^M^; /= '"^- (37) 
|fc'||9'l 

The entries of the matrix 'E^ (W, W , x') are given in appendix Apphcation of Eqs. ( p6| ) and ( pTj ) leads to a system 
of coupled integral equations for the multipole amplitudes, 

ReM'i{W)^Mi'^°^^{W) + -V ( dW ^JCii,{W,W')luiM\{W'). (38) 

^ -ir 

The integral kernels 

1 

lC{i,{W,W')^ j dxT>i{x)li{W)E\W,W' ,x')Yr^{W')Gv{x') (39) 
-1 

couple the multipoles with each other. Analytic expressions of these kernels have been given in Refs. §,0 for the 
S and P waves. These kernels are regular kinematical functions except for the diagonal kernels ]C\i, which contain 
a term cx 1/(VF — W). Our further analysis is based on Eqs. (^8|), which may be solved by several methods. In j|] 
conformal mapping techniques have been used. The work of Schwela and Weizel Q follows the method of Omnes 
pTf . Lebedev uses Pade approximants. Since our aim is to analyse new experimental data and not to give mere 
predictions, we essentially follow the work of leading to a natural parametrization of the multipoles. Before we 
turn to this point in the next section, we have to make some remarks. 

The dispersion relations (|33| ) are valid for the isospin components with / = ±,0. Since the unitarity condition (|2^ ) 
will be imposed on the solutions for the multipoles, we have to work with amplitudes of definite isospin i or |. This 
leads to a mixing of the components with / = + and / = — . The corresponding linear combinations of the integral 
kernels are obtained as follows. If the kernel connecting the multipoles M.\, and Ai\ (I = 0, ^ | in the following) is 



denoted by /C/;{ , Eqs. (18) lead to the relations 



3 1 
2 • 3 
1 ..jr 1 



/' = -: /Cf/ =-(4/C+ +c^/C„,), (40) 
^' = : IC^, = CqIC"i, , 

with c^"^ = c^^ — Cq = 1, cl^ = — 1, cl"^ = 2 and all other coefficients vanishing. 

In this work we restrict ourselves to the energy region up to « 500 MeV. The first reason for this are possible 
problems with the convergence of the partial wave expansion which might occur in the unphysical region, when the 
path of the dispersion integrals leads outside of the Lehmann ellipses |13| . Furthermore, in this energy region it is 
reasonable to restrict the solution of the integral equations to S, P and eventually D waves. At energies near the 
second resonance region, higher partial waves have to be taken into account which will considerably complicate the 
procedure. Finally our method provides the most reliable results as long as the complex phases of the partial waves 
are known. This is the case below two-pion threshold. Nevertheless it can be assumed that the effects of inelasticities 
remain small also at somewhat higher energies. As a consequence of these restrictions, our results will be especially 
useful for the study of low energy amplitudes and the electromagnetic excitation of the A(1232). 

IV. SOLUTION OF THE INTEGRAL EQUATIONS 

On condition that the complex phases of the solutions are known, Omnes lO] has proposed to solve the integral 



equations (38) as follows. Introducing the functions 

hl{W)^e^^'^'^Uin(pl{W), (41) 
with the complex phase (f)\ [W) of the multipole M.\, one obtains 

ImMiiW) = h\*{W)M\ {W). (42) 
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Next one defines the amplitudes 



MjiW) ^ri{W)Mj (43) 
by taking out kinematical factors for the S and P wave amplitudes, 

W 



D,Df{W-mN)' 



'^1- ^ T^, ^ ' (44) 



ri+{W) 



k I {W-mN)D, 
W 

D,Df{W-mN) \ k\\q 



with Dij — y/Ei_f + m^v. The coupled system of integral equations is now cast into the form 

hl*{W')Ml{W')dW' 



Mf(W-)=MfP°l^(Ty) + i J 



W'-W-ie 
^thr 



OO 



+-E / K''iW,W')hi;*iW')Mi:{W')dW', (45) 
' '^thr 

where the singular term has now been treated individually and subtracted from the integral kernel, 

In order to treat Eqs. (^5|), some further simplifications and assumptions are necessary. First of all, the infinite sum 
over /' has to be cut off. It can be shown that for multipoles with small complex phases in the relevant energy region, 
the solutions of the integral equations are essentially given by the inhomogeneities, i. e. the pole terms. Therefore we 
solve Eqs. (^ ) only for the S, P and £>i3 waves, and represent the higher partial waves by the inhomogeneities. A 
study of the integral kernels shows that many of them can be neglected. This leads to further simplifications. 

To solve Eqs. (45) the complex phases 4)\ {W) of the multipoles have to be known over the whole range of integration. 



From the Fermi- Watson theorem ( |28D we know that these phases are equal to the pion-nucleon scattering phase shifts 
below two-pion threshold. Above the inelastic threshold, 4)\ {W) is no longer uniquely determined by unitarity. 
However it has been shown that the following ansatz for (f)\ {W) can be motivated from the unitarity relation of the 
inelastic amplitudes 

(W) = arctan , ' ^ , ^ ./ \ ' . (47) 

Below two-pion threshold the ansatz (|7|) leads to the vrA^ scattering phase shifts 5j. Above this threshold the 
inelasticity parameters rjf as determined from pion-nucleon scattering enter into the formula. The ansatz (^ ) is 
justified if both 7iV — s- nirN and TriV — s- wkN are dominated by the same resonant intermediate state or if both 
matrix elements are real (e.g. if they are both described by one particle exchange) or imaginary (e.g. in the case 
of purely diffractive processes). In the case of the P13 waves the ansatz (47) would lead to an unphysical resonant 
behaviour. Therefore in this case we choose the alternative ansatz 



<b'(W) - arctan ( vHW) ^ra2S{ jW) 



(48) 



which has also been derived in Ref. |^ . 

At higher energies, where information from pion-nucleon scattering is no longer available, certain assumptions have 
to be made in order to interpolate the phases to infinity ||^ . As an alternative we cut off the integrals of Eq. ( p5| ) at 
W = 2 GeV and represent the high energy tails of the dispersion integrals by t-channel exchange of p- and cj-mesons, 
leading to the contributions: 
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my zm^r t — my my ztun t — rriy 



ty • t - my 



where V stands for w or p respectively, Ay is taken from the radiative decay of the vector mesons and gy are the 
vector and tensor couphngs at the VNN vertex. The integral equations now take the form 

^thr 



A 

J K'iW,W')hi:*{W')Mi:iW')dW' + Mi^^{W), 



with a cut-off at A = 2 GeV. The high energy parametrization by t-channel exchange leads to a further inhomogeneity 

of the integral equations, m['^(VF), to be constructed according to the rules given above. In comparison with jl^] 
this leads to a slight modification of the solutions. As in we first solve the characteristic equations 



^ h{*{W')Ml{W')dW' 



W ~W -ie 
^thr 

+m\^{W), (51) 
which contain only the singular part of the integral equations. The general solutions are 

Mf {W) = Mf'P°^^(PF) + m{^{W) 

Mf^°"^(W^) f h\{W'){M\'^'''^''{W') + M\^{W'))dW' 



thr 



m\^'''^{W'){W' -W -ie) 



with 



+P(VF)m[;J°™(1¥), (52) 



../,hom.„.. I W f ct>iiW')dW' I 

= V J W'iW'-W-^e) (''^ 



^thr 



a solution of the homogeneous equation 



^hr 

and P{W) an arbitrary polynomial. In our calculations for this polynomial only the simplest case, i.e. a real constant, 
will occur. 

If the self-coupling of Mf (W) is fully taken into account, Eq. ( pO[ ) has to be transformed into an equation of the 
Fredholm type, 

A 

Mj{W)=Nl{W)+J2- f Fjl'{W,W')hl'*{W')Mi:iW')dW', (55) 



^thr 
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with 



N^iW) = Mf'P°^°(T4^) + m\^{W) (56) 



^thr 

and the Fredholm kernels 



W^') ^ K\{:{W,W')hl'*{W') + 



A 



/if (T4^")Kf,{' (PF", W')h\* {W')dW" ^^^^ 

- Mf;^°™(W^")(Ty"-M/-ze) 
tnr 

Eq. (HI) can now be solved by standard numerical methods. We mention that also the solution of Eq. (^) is not 
unique. A solution of the corresponding homogeneous equation can be added in the same way as in the case of Eq. ( |5l| . 
For the mathematical details of solving the homogeneous equation see Ref. . 

Numerically we found that the self-couplings of the multipoles are given to good approximation by the singular 

(-) (-) 

parts of the integral equations, except for the P33 multipoles Mj^^ and E]^ . Therefore, the Fredholm equations have 
to be solved for these P-waves, while in all other cases it is quite sufficient to solve the characteristic equations. If 
the couplings to other multipoles are taken into account, they enter into the equations as further inhomogeneities. If 
we neglect all integral kernels giving only small contributions, we obtain the following procedure to solve the integral 
equations for the S-, P- and £'i3-waves: 

(3) . . . 

1. The amplitude M-^^ is practically independent of the other multipoles. The corresponding Fredholm equation 
has to be solved. 

(-) (-) 

2. Having solved the equation for M{^\ we calculate its contribution to E]^ entering into the corresponding 

equation as a further inhomogeneity. 

(-) (-) (-) (-) 

3. The equations for Ml_ and Ml_ are solved taking account of the influence of Mj^^ and E]^ . 

(-) (-) 

4. For £'0+ and i?o^ also the contributions of the Mi_ components are important. 

5. In the isospin channel, only the contribution of M^} to E^_^ has to be taken into account. 

6. For all other P waves as well as the Dia multipoles we only have to solve the characteristic equations. 



According to Eq. (|5^) we can add to each solution of the inhomogeneous equation a solution of the homogeneous 

equation multiplied by an arbitrary polynomial in W. Following the arguments of we restrict ourselves to the 

(0) (-) fo) (-) (-) {-) 
addition of homogeneous solutions to the multipoles E^^, Ei^^ , M^J, M^t , i?^^ and M^^ multiplied by a real 

constant. These six constants are determined by a fit to the data. Further parameters of our fit are the coupling 

constants of the vector mesons. In Q it is argued that vector meson exchange should not be included because of 

double counting. This argument is correct because the authors evaluate the occuring integrals to infinity. From a 

study of unitarity diagrams it can be seen that the integrals in fact contain the vector meson contributions. In our 

case, however, the integrals run only to the cut off' A, and only a fraction of vector meson exchange will be contained 

in the integrals. The remaining part is assumed to be the dominant contribution at the higher energies, which we 

represent by vector meson exchange as described with four free VNN coupling constants fitted to the data. Thus we 

end up with ten free parameters. 

In solving the integral equations we have to know the scattering phase shifts and inelasticity parameters of pion 

nucleon scattering. They are taken from the SAID code (solution SM95 jl^). 
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V. THE FIT TO THE DATA 



We have determined the ten free parameters of our approach by a fit to selected photoproduction data for 160 MeV 
< < 420 MeV (see Tab. |). Most interesting for our analysis are the new data from MAMI. Differential cross 
sections of tt" production off the proton near threshol d |l5| , and differential cross sections as well as beam asymmetries 
E for both tt"*" and tt" production off the proton | |T^ , |l'<| have been measured with high precision. Since there are no 
new low energy data for Tr~^ we include those of the data compilation ||l^ . We have also included the most recent data 
from ELSA, differential cross sections for 7r+ production |l9| and the target asymmetries T for 7r+ and tt*^ production 
off the proton |2^. To get the full isospin decomposition we also had to include data for tt^ production, for which we 
took the differential cross sections given in Refs. pT|-p5t. Altogether we have included 1751 data points. The overall 
per data point of 2.8. We note that we have not renormalized the different data sets in our fit. 

In Figs. |i|-p^ we show our results for the observables. In general we obtain a good description of the photoproduction 
data. The differential cross sections and beam asymmetries for 7r+ production are in perfect agreement with the Mainz 
experiment except for the differential cross sections around 330 MeV where we find a slight discrepancy (see Fig. |^). 
In our local fit (see next section) this discrepancy could not be removed which means that there is a problem either 
with our parametrization or with the experiment. We stress that also at these energies the beam asymmetries, which 
are less sensitive to systematical errors, are well described by our fit (see Fig. ||). From the experimental point of 
view one main difficulty are measurements concerning 7r+ production at very low energies. Because of the short life 
time of charged pions many of them decay before detection. As a consequence the existing data have large error bars. 
In Fig. |l| we show the result of the fit compared to the world data. We find good agreement within the statistical 
uncertainties. The differential cross sections and target asymmetries from Bonn which we included in our fit (see 
Figs. 1^ and ||) have larger errors than the new Mainz data. In the case of the target asymmetries our result seems to 
be consistent with the data if the statistical uncertainty is taken into account, except for very low energies where the 
experimental values are systematically larger than our curves. However, it can be shown that such large values of T 
would require much larger imaginary parts of the partial waves than allowed by unitarity considerations. 

In tt" production on the proton we find a somewhat different situation. In this case the results of recent high 
precision experiments at low energies are available ||l^. However, only the differential cross sections for E-y > 160 
MeV were included in our fit. The description is reasonably good also for the lower energies (see Fig. ||). At higher 
energies we again find a very nice description of the beam asymmetries (see Fig. |^). Also the differential cross sections 
are well described except for the energies around 330 MeV (see Figs. M and 0). The target asymmetries are only in 
semiquantitative agreement (see Fig. [lO|) . 

In pion production on the neutron the experimental situation is less satisfactory than for the proton. The most 
precise experiment which has been performed at TRIUMF by studying the inverse reaction tt~p — > 771 psf , is restricted 
to energies below the A resonance. In Fig. |ll| it can be seen that we fit this experiment very well. However, to take 
account of a larger energy range, we had to include older data which are reasonably well described within the larger 
statistical uncertainties (see Fig. |l2|). 

The parameters determined by our fit have no direct physical interpretation. Given the integral equations, the 
coefficients of the homogeneous solutions have to be determined by the boundary conditions determined by the physics. 
Concerning the coupling constants of the vector mesons, we expect them to be of the same order of magnitude as 
the corresponding values of other analyses. In Tab. || we compare our results with other references. The existing 
differences should be attributed to the fact part of the vector meson strength is already contained in the dispersive 
integral up to 2 GeV. 



VI. RESULTS, PREDICTIONS AND DISCUSSION 



A. S and P wave multipoles 

In Figs. [T^-p^ we show the isospin components of the S- and P-wave multipoles with / = 0, The solid 

lines show the results for the real and imaginary parts of the "global fit" obtained with the complete data set. The 
error bands (dashed lines) result from the errors in the parameters as determined from the covariance matrix after 
normalizing to 1. Because of the relatively small number of parameters, the errors of the global fit are very small 
except for the components of i?o+ and Mi_ with / = 0, which reflects the uncertainty of the neutron data. In order to 
visualize the fluctuations of the data, we also analysed them in bins of 20 MeV with the same parametrization ( "local 
fit"). Except for 7r° production on the proton, the data basis for energies E^ < 210 MeV is quite unsatisfactory. As 
a result we find the largest error bars for all amplitudes in this region. Figs. pT|-po| show the decomposition of the 



isospin i amplitudes into proton and neutron part according to Eqs. (19). In this case, the corresponding local fit 
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was obtained by first fitting the pM^^) and M'^^ components to the proton data, and then determining the „M*^^) 
amplitudes by a fit to the neutron data. In the case of the proton, the local fit is in general agreement with the result 
of the global fit, while we find strong fluctuations and large error bars for the neutron amplitudes. It is obvious that 
more and better neutron data are needed to get full information on the amplitudes £^0+ and A/i_ in the isospin ^ 
channel, which are of special interest for studies of the associated resonances 5ii and Pn. 

Given the results for the S and P waves from our fit, we now turn to a detailed study of the threshold amplitudes 
for both charged and neutral pions and to the electromagnetic excitation of the A(1232). 

B. Low energy amplitudes 

As pointed out before, our data basis does not include the threshold region. The main reason for this is that we have 
to expect isospin symmetry breaking in that region due to the pion mass splitting, while our dispersion calculation 
is based on the isospin symmetry of the pion-nucleon phase shifts. Therefore, the threshold S wave amplitudes given 
below are a prediction, in the sense that the cross sections above E-^ — _L60 MeV determine the threshold values by 



continuation of the analytic functions constructed in dispersion theory JM] . Table III compares our threshold values 
for the charged pion channels to the "classical" low energy theorem |p7| (LET), chiral perturbation theory ^8 29 
(ChPT) and experiment. Note that ChPT contains lowest order loop corrections, while LET is based on tree graphs 
only. The agreement between our calculation and the experiment is excellent. 

It should be noted that a larger value for i?o+(7r^p), namely — 34.7 ± 1.0 (here and in the following in units of 
10^'^/to^), was reported |]30| . However, this value derived from angular distributions was recently withdrawn 
Incidentally, that (wrong) value would have resulted in a difference of ai — 03 = 0.275/m7r for the pion-nuclcon 



scattering lengths, while the (correct) value of Tab. Ill gives a smaller value of (0.253 ± 0.003)/m7r. However, the 
larger value would be in much better agreement with the analysis of pion-nucleon scattering by the Karlsruhe group 
101 and recent results from pionic atoms which lead to a value of (0.276 ± 0.013)/m7r |^^. Since this cross relation 
between hadronic and electromagnetic physics is only based on the well-known Panofsky ratio, a{'K~p ir'^n) / 
a(TT~p — *■ "fn), even the small deviation in oi — 03 poses a problem. Fig. |2^ compares our results to the angular 
distributions near threshold measured at TRIUMF |Q. With the exception of the data at 165 MeV we find excellent 
agreement. 

Neutral pion photoproduction near threshold has been a topic of many experimental and theoretical investigations 
in the 90's. The classical LET predicted a power series expansion of the amplitude in /i = m-^/mN ~ 1/7, leading 
to a threshold value of about —2.4 plus presumably small corrections in fi^. The experimental findings at Saclay 
p3[ , Mainz and now at Saskatoon |^ are a considerably different amplitude of about half that value. The 
classical LET was based on Lorentz, gauge and chiral invariance It further assumed that the amplitudes are 

analytical functions, though In fi singularities due to pion loops had been occasionally discussed in the literature. 
The first consistent calculation of pion photoproduction to one-loop order in ChPT clearly showed the breakdown of 
analyticity. The result is the expansion |^ 

I _0\ ^9ltN ( 2 3 + 'tp 2 "^W , r 3A lro\ 

) = - ^ - ^ 1671 + ] j ' ^ ^ 

with the first two terms on the rhs given by the old LET, and the third term a correction due to pion loops. With 
the numerical values for the anomalous moment Kp, the nucleon mass m^r, and the pion decay constant = 93 
MeV, this correction even changes the sign of the threshold amplitude. Various predictions are compared with the 



experiment in Tab. HI. The value given by ChPT was recently evaluated in the heavy baryon formalism up to order 
|29| . Our value from dispersion theory is obtained by breaking isospin symmetry explicitly by the physical masses 
of the pions and by assuming that the imaginary parts of the amplitudes start at charged pion threshold and are 
given by the results of our global fit, i.e. by iteration of Eq. (|5^) under the stated assumptions. The agreement with 
experiment is astounding. Once more the data at the higher energies "predict" the threshold values with a good 
accuracy. This result is particularly surprising, because the real part of _Bo+ obtains very large distributions from the 
imaginary parts of the higher multipoles via the dispersion integrals. Altogether these contributions nearly cancel the 
large contribution of the nucleon pole, which corresponds to the result of pseudoscalar coupling, leading to a total 
threshold value 

Re£;5^'^(p7r°) = -7.63 + 4.15 - 0.41 + 2.32 + 0.29 + 0.07 = -1.22, (59) 
ReE*^'^(n7r°) = -5.23 + 4.15 - 0.41 + 3.68 - 0.93 - 0.05 = 1.19, (60) 
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where the individual contributions on the rhs are, in that order, the pole term, Mi+, £^0+, Mi_ and higher 
multipoles. We add that the existing uncertainties for the neutron value of Mi_ give rise to large effects. Fig. ^ 
compares our results to the measured angular distributions and predicts values for the three single polarization 
observables. The dotted curves are obtained by arbitrarily reducing Mi_ by 50 %, which demonstrates the strong 
sensitivity of E and P with regard to that multipole. Finally, Fig. shows the energy dependence of i5o+ in the 
threshold region. The existence of a Wigner cusp in Rei?o+(p7i'°) at charged pion threshold is due to the infinite 
derivative of Imi?o+ at that threshold, and a reflection of the coupling to the much stronger 7r+ channel. We also 
calculated the threshold values for the three P waves and find agreement with the prediction of ChPT and experimental 
data within 5-10 %. For a detailed discussion see Ref. | p6| |. 



C. The ratio E1+/M1+ in the A region 



The search for a deformation of the "elementary" particles is a longstanding issue. Such a deformation is evidence 
for a strong tensor force between the constituents, originating in the case of the nucleon from the residual force of 
gluon exchange between the quarks. Depending on one's favourite model, such effects can be described by d-state 
admixture in the quark wave function, tensor correlations between the pion cloud and the quark bag, or by exchange 
currents accompanying the exchange of mesons between the quarks. Unfortunately, it would require a target with a 
spin of at least 3/2 (e.g. A matter) to observe a static deformation. The only realistic alternative is to measure the 
transition quadrupole moment between the nucleon and the A, i.e. the amplitude While the existence of such 

a transition cannot be regarded as direct proof of a deformation, it is sensitive to model parameters responsible for 
possible deformations of the hadrons. 

The experimental quantity of interest is the ratio Rem = in the region of the A. The quark model with 

SU(6) symmetry predicts Rem = 0. Depending on the size of the hyperfine interaction and the bag radius, broken 



SU(6) leads to —2% < Rem < |3^j3^. Similar results can be obtained in the cloudy bag model |39|j40|], larger 
effects with Skyrmions (—5% < Rem < —2%). Exchange currents yield values of about —3.5%. The first 
lattice QCD results are (+3± 9) % The analysis of the data changes by time and method, with a PDG average 
of (-1.5 ±0.4) % 0. 

The comparison of our dispersion analysis with the analysis of the Mainz data by Beck et al. shows nice 
agreement with the multipole Mi^ but systematic differences for the quadrupole transition Ei^ at energies above 
the resonance. However, the ratio Rem at resonance agrees quite well in both analyses (see Fig. p3). Unfortunately, 
there appears to be a discrepancy of up to 15 % between the differential cross sections measured at Mainz (upon 
which our analysis is largely based) and the LEGS data ^] at Brookhaven. As far as ratios of cross sections are 
concerned, the data seem to be in general agreement. As an example, the ratio of the LEGS cross sections for 
photons with polarization parallel and orthogonal to the scattering plane , is nicely described by our analysis (see 
Fig. p5| ), though these data have not been included in our fit. In view of the importance of the quantity Rem for 
our understanding of the internal structure of the nucleon, the existing discrepancies between the data deserve further 
studies. Based on "our" data selection we obtain a value Rem = (—2.5 ±0.1) % for the global fit ((—2.33 ± 0.17) % 
for the local fit), which is close to the value derived by Beck et al. p9| . 

Another point of interest is the separation of resonance and background contributions. The problem in this quest 
arises, because the leading tree graphs, the Born terms and the A excitation, do not fulfill the requirements of unitarity 
(Fermi- Watson theorem). In physical terms, the pion produced by the Born graph feels the presence of the resonance 
and is reabsorbed by the nucleon to form a A, which altogether leads to the required phase relation of Eq. (|2^). 
Together with the "pure" A, these diagrams constitute the "dressed" A. A technique to achieve this result is the 
method of Olsson [ p0[ , 

Afi+ = |A/f+|e'*« + |Afi7|e*''^<=^e** = \Mi+\e'^^+ , (61) 
which combines the two contributions, unitarizcd individually, by a rotation of the resonance contribution in the 



complex plane. The other technique is the K- matrix method 51 



K{W) = A{W)/iMr,s -W)+ B{W), (62) 



T{W) = K{W)cosS{W)e'^'-^\ (63) 

where B(W) is a slowly varying background, and S{W ~ Mrcs) = tt/2. It is obvious that different approaches will 
give different answers to the question which part of ratio Rem is due to the "pure" resonance. 
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Due to the analyticity of the functions involved, dispersion theory is particularly apt to study the amplitudes in 
the complex energy plane. Therefore, we have analyzed the amplitude by comparing the "speed" |5^ ] 



SP [M,+ iW)] 



dMi+(W) 



dW 



(64) 



with the speed of an ideal resonance due to a pole in the complex plane at Wr — Mi^ — iV ji/2. Such a resonance has 
the amplitude 

'g'''"''- A/.-t'-"r„/2 ■ 

with re** determining the complex residue at resonance, W = Wr. 

Due to the differentiation involved, the effect of a slowly varying background is much reduced, and the "experimen- 



tal" speed agrees nicely with the "ideal" speed, as can be seen in Fig. 26 (left part). The same is true for the speed 
vector, the complex number dAIi^/dW , which has been plotted on the rhs of Fig. Gn. The experimental values given 
by the crosses in this figure agree quite well with the squares of the ideal resonance, except for the smallest energies 
close to the origin. We note that as function of energy W the values run counter-clockwise, starting near the origin at 
threshold (the "stem" of the apple-shaped structure), passing through resonance at maximum speed (the "blossom" 
of the apple), and again approaching the origin for ^ oo. The fact that the "apples" for Mi+ (top) and Ei+ 
(bottom) are essentially oriented in opposite directions is related to the negative value of Re m- The phases of the 
complex amplitudes at resonance are defined by the angles and respectively (see Fig. |26|) . Table |^ shows 
our results for the pole position and the complex residue at the pole for both the magnetic dipole and the electric 
quadrupole transitions. The complex ratio of the residues is given by 

Rem = ''''^If^^l = -0-035 - 0.046i. (66) 

The values for the pole position are in perfect agreement with the analysis of pion-nuclcon scattering. Essentially 
based on analyses of Hohler et al., the PDG ||| hsts Mr = (1210 ± 1) MeV and Tr = (100 ± 2) MeV. The complex 
ratio Rem ^ pole, W — Mr — irR/2, is rather stable in comparison with the value Rem on the real axis, at 

W = Mj-cs ~ 1232 MeV. As an example, Sato et al. have studied two versions of a dynamical model leading to a 
difference in Rem '-'^ ^0 % while Rem differs by as much as 300 %! It is also worth noting that studies of hyperon 
resonance radiative decay in the framework of ChPT have predicted values for Rem compatible with our results 
within the error bars. Finally, Fig. ^ shows the resonance-background separation of the real and imaginary parts for 
the two A multipoles. Indeed, the speed-plot technique leads to a slowly varying background except for the threshold 
region where the background has to cancel the "ideal" resonance shape to give the (vanishing) physical amplitude. 
While the experimental amplitude particularly its real part, is quite different from a typical resonance shape, 

the amplitude after background subtraction is not too far from an ideal resonance. It is in fact rather astounding 
that Mr can be predicted to four digits and Tr to two digits from the £'1+ amplitude, which is about half resonance 
and half background (see Fig. |27| ). 

Instead of the multipoles the electromagnetic resonance excitation is very often discussed in terms of helicity 
amplitudes. In Tab. ^ we give the helicity matrix elements Ai and for the transition jN A which are derived 
from our results for the multipoles via the standard definitions. In that table also the ratio Rem is given which is 
equivalent to the result of the application of the K-matrix method. 



D. Total cross sections 



As a result of our fit we also obtain predictions for the total cross sections. These are shown in Fig. |2^ for the 
four physical reactions, together with the contributions of the leading S- and P-wave multipoles. In terms of the 
multipoles, the total cross sections are 

I h I 

^ 2^i-i ^(Z + 1)2 [(/ + 2) (I p + I Af,+i,_ p) + I (I M,+ p + I p)] . (67) 

' * ' /=o 

with the isospin combinations of the multipoles according to Eqs. ([l^-[l7|). 
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A topic of much current interest is the spin structure of the total absorption cross section, which will be measured 
with circularly polarized photons and nucleon target polarization in the direction of the photon momentum |^ , ^ . 
Depending on their relative direction, the spins of photon and nucleon add to spin 3/2 or subtract to spin 1/2, which 
defines the helicity structure of the total cross sections, era and ai, with a^^f. ~ {as + cri )/2. The difference of these 
two qualities enters into the Gerasimov-Drell-Hearn (GDH) sum rule, 

CO 

which relates the helicity structure of the total absorption cross section of the nucleon with its ground state properties. 
Our prediction for the contributions of single pion production to as — ai and the GDH integral as function of its 
upper limit is shown in Fig. |2^. 



VII. SUMMARY AND OUTLOOK 



Dispersion relations at fixed t provide a useful framework to analyze pion photoproduction between threshold and 
an excitation energy of about 500 MeV. After multipole projection we obtain a system of coupled integral equations 
for the multipole amplitudes is whose solutions depend on 10 free parameters, which are fitted to selected experimental 
data in the range of 160 MeV < < 420 MeV. On the basis of these amplitudes, we obtain predictions for the 
threshold behaviour of pion photoproduction, including the cusp effect for neutral pion production. The results for S 
and P waves agree within about 5 - 10 % with the predictions of ChPT and a direct analysis of the threshold data. 

In the region of the A (1232), we find a ratio of the electric quadrupole to the magnetic dipole amplitude of 
Rem = (—2.5 ± 0.1)% at W = M^^s = 1232 MeV. The A resonance pole in the complex plane can be determined 
from both amplitudes, in excellent agreement with the results of pion-nucleon scattering. The ratio of the complex 
residues at the pole is R^m ~ (—3.5 — 4.6i)%. While Rem is very model dependent, the value at the resonance pole 
is relatively stable. 

Since our "global fit" contains only relatively few parameters, it cannot follow all the fluctuations of the data in the 
energy and angular distributions. Therefore, the overall is only 2.8 per data point. We are convinced that many of 
the fluctuations of the data are due to systematical errors and, therefore, should not be reproduced by a theoretical 
analysis. In order to obtain an estimate of such systematical errors, we also show the result of a "local fit" by fitting 
the free parameters to, e.g., the data in 20 MeV bins individually. 

As further tests of the remaining model dependence of our procedure, we are presently studying some modifications 
of the parametrization. First, the high i-channel contributions to the dispersion integrals will be replaced by a Regge 
model. Second, all isospin components of the S, P and D waves will be treated on the same footing by allowing 
for the addition of homogeneous solutions to all of them. The explicit fit of all D wave contributions will give the 
possibility to study the influence of the higher partial waves. 

In conclusion we have obtained a good representation of the data in the range up to about 500 MeV. The biggest 
experimental uncertainties are shown to be due to the neutron data, in particular in the case of the Eqj^ and Mi_ 
amplitudes. New experiments to determine the small amplitudes by polarization observables are urgently needed, 
particularly in the case of the neutron and the Roper multipoles. 

APPENDIX A: MATRICES CONNECTING THE AMPLITUDES 

The matrices connecting the various amplitudes defined in the text are: 



H = diag 



I ^(g/+mw)(g,+mw)(W-™iv) \ 
-y/(-E.+mjv) (vy-mjv)|fc| 



^(£;j+mjv)(g.+mjv7(VF-mjv)|fc| 
ygj+mw (VK-rrajv)|fc|^ 



(AI) 
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/ 1 

-1 



V 







W — TTlN 

-(W + mjv) 



2{W-mN) 
t 



W — mjv + 



2{W-m_w) 
t 



4W^ 



2W 

2W 2W {C-%3 

\ 2W 2W W 



W—rriN 



-2W 

W+niN 



with 



{C-^)ii=2W{W + mN), 
{C-^)i2 = -2W{W-mN), 
t-ml 



(C-i)i3 = 2mNW 
(C-i)i4 = 2mNW 



W — 
t — 

W + tjin' 



{C-^)SA = 2W -TUN 



2{W-mN) 
2{W + mN))' 



Gi{x) 



Pl^.ix) PU{x) lPl+,(x) (l + l)Pl_,{x) 



(; + l)P/(x) 
-Pl'ix) -P['{x) Pl'ix) 



PlUix) Pl'_,{x) 



IPlix) 



D((x) 



UP' -P'-^ ^(Pl+l-Pl-l) ^,{Pl-Pl-2)} 
' {Pi -Pl+^ 2kliP+^-P-^) 0} 

0} 



Pz-l) 

2(7+1)1.-'' -.1+1 2i+TV-'H-i Pl-l) 
\ i{-P, Pi-, ^(P,_i-P,+i) 



with the Legendre polynomials P; as functions of x. 
The entries of the matrix E are: 



El, 



^12 



^13 



Eli 



-^21 



^22 



-2^23 



1 



+ 



W'-W W'^-u 



+ 



W — niN 



W-W + 



W — mjv 



W' + W W^-u 

l-e^ t-ml 

2 {W - mN){W - ruNY 

1 — £^ t — m-l 

2 {W +mN){W -mN)' 

1 t-ml \ 

1 I W — W -\ — 

W + W W'^ -u\ W + mNj 



1 



w 



(w' + w+ '-"^^ ] 

W'^-u V W + mNj 



t-mi 



2 {W -mN){W + mN)' 



15 



^24 = 



I 



t-mi 



-E'32 



-^33 



^34 



^42 



-2^43 



^44 



2 (VF'+TOjv)(M^ + mAr)' 



1 

+ 



W'-W W'^-u 
1 



2s' 



2e' 



+ 



+ 



W'-W w^- 



W-W' - 2mjv 



W — mjv J ' 



-W^' - + 2mjv - 



W + mjv 



-W^' - - 2mjv - 
W -W' + 2mN 



t-nil 
W'-ruN 
t-ml 



W' + mN J ' 



(A6) 



APPENDIX B: RESPONSE FUNCTIONS 

The response functions discussed in this contribution may be expressed by the CGLN amphtudes as follows: 

iir = I jFi |2 + I jF^ P +lsm^e{\ J'al^ + l T4 I') 

-Re{2 cos e TlTi - sin^ Q {TIT4 + + cos 6» TIT^)}-, 

i?^' = sin 6* \wL{TlT-i - T^T^ + cos6» {TIT4 - T^T-i) - sin^ 6* TIT^S, 

i?7 = - sin 61 Im{2jr*jr2 + TlT-i - T^^T^ + cos Q {TIT 4 - T^Tz) - sin^ 6 T^Ta}, 

Rtt = sin^ e (i(| ^3 I' + I Ti n + Re{^r^4 + T^T^ + cosO T;Ta}). 



APPENDIX: ACKNOWLEDGMENTS 



We would like to thank Prof. G. Hohlcr for very fruitful discussions and the members of the A2 collaboration at 
Mainz for providing us with their preliminary data, in particular J. Arends, R. Bock, F. Hartcr, H.-P. Krahn and 
H. Stroher as well as G. Anton from ELSA at Bonn. This work was supported by the Deutsche Forschungsgemeinschaft 
(SEE 201). 



[1] G. Hohler, Pion-Nucleon Scattering, Landoldt-Bornstein, Vol.I/9b2, cd. H. Schopper, Springer (1983). 
[2] I. S. Barker, A. Donnachie and J. K. Storrow, Nucl. Phys. B 95 (1975) 347. 
[3] D. Schwela and R. Weizel, Z. Physik 221 (1969) 71. 



16 



G. F. Chew ot al., Phys. Rev. 106 (1957) 1345. 

J. D. Bjorken and S. D. Drell, Rolativistic Quantum Mechanics, McGraw-Hill, New York (1964). 

K. M. Watson, Phys. Rev. 182 (1954) 228. 

D. Drechsel and L. Tiator, J. Phys. G 18 (1992) 449. 

J. S. Ball, Phys. Rev. 124 (1961) 2014. 

F. A. Boronds, A. Donnachie and D. L. Weaver, Nucl. Phys. B 4 (1967) 1. 
O. Hanstein, Ph.D. thesis, Mainz (1996). 
R. Omnes, Nuovo Cim. 8 (1958) 316. 

A. I. Lebedev und B. V. Mangaaeyev in A. A. Komar (ed.), Photoproduction of Pions on Nucleons and Nuclei, Nova 
Science Publishers, New York and Budapest (1989). 

A. Donnachie in E. U. S. Burhop (cd.). High Energy Physics, Vol. V, Academic Press, New York (1972). 
R. A. Arndt, I. I. Strakovsky and R. L. Workman, Phys. Rev. C 52 (1995) 2120. 

M. Fuchs et al., Phys. Lett. B 368 (1996) 20, 
M. Fuchs, Ph.D. thesis, GielSen (1996). 

H. -P. Krahn, Ph.D. thesis, Mainz (1996), 

R. Beck, H.-P. Krahn et al., Phys. Rev. Lett. 78 (1997) 606. 
F. Haerter, Ph.D. thesis, Mainz (1996). 

D. Menze, W. Pfeil und R. Wilcke, Compilation of Pion Photoproduction Data, Bonn (1977). 

K. Biichlcr et al., Nucl. Phys. A 570 (1994) 580. 
H. Dutz, Ph.D. thesis, Bonn (1993), 
D. Kramer, Ph.D. thesis, Borm (1993), 

B. Zucht, Ph.D. thesis, Bonn (1995), 

H. Dutz et al., Nucl. Phys. A 601 (1996) 319.. 
M. Boncvcntano et al., Nuovo Cim. 4 (1958) 323. 

F. Carbonara et al., Nuovo Cim. 13 A (1973) 59. 
P. Benz et al., Nucl. Phys. B 65 (1973) 158. 

T. Fuji et al., Nucl. Phys. B 120 (1977) 395. 
A. Bagheri et al., Phys. Rev. C 38 (1988) 875. 

O. Hanstein, D. Drechsel and L. Tiator, Phys. Lett. B 399 (1997) 13. 
P. de Baenst, Nucl. Phys. B 24 (1970) 633. 

V. Bernard, N. Kaiser and U.-G. Meil3ner, Phys. Lett. B 378 (1996) 337. 
V. Bernard, N. Kaiser and U.-G. Meifiner, Phys. Lett. B 383 (1996) 116. 
K. Liu, Ph.D. thesis. University of Kentucky (1994). 
M. A. Kovash ot al., ttN Newsletter 12 (1997) 51. 
D. Sigg et al., Phys. Rev. Lett. 75 (1995) 3245, 

D. Sigg et al, Nucl. Phys. A 609 (1996) 269, 

M. Janousch et al., Proc. of the 14*'' Int. Conf. on Particles and Nuclei, Williamsburg 1996, eds. C.E. Carlson, J.J. 
Domingo, World Scientific, Singapore 1997, p. 372. 

E. Mazzucato et al., Phys. Rev. Lett. 57 (1986) 3144. 
R. Beck et al., Phys. Rev. Lett. 65 (1990) 1841. 

J. C. Bergstrom, Phys. Rev. C 52 (1995) 1986, 

J. C. Bergstrom et al., Phys. Rev. C 53 (1996) 1052, 

J. C. Bergstrom, R. Igaraslii, and .J. M. Vogt, Phys. Rev. C 55 (1997) 2016. 
V. Bernard, J. Gasser and U.-G. Meifiner, Phys. Lett B 268 (1991) 291. 
R. Koniuk and N. Isgur, Phys. Rev. D 21 (1980) 1868. 
D. Drechsel and M. M. Giannini, Phys. Lett. 143 B (1984) 329. 

G. Kalbermann and J. J. Eisenberg, Phys. Rev. D 28 (1983) 71, and 
Phys. Rev. D 29 (1984) 517. 

K. Bermuth, D. Drechsel, L. Tiator and ,J. B. Seaborn, Phys. Rev. D 37 (1988) 89. 

A. Wirzba and W. Weise, Phys. Lett. B 188 (1987) 6. 

A. Buchmann, E. Hernandez and A. Faessler, Phys. Rev. C 55 (1997) 448. 

D. B. Leinweber, et al., Proc. Int. Conf. "Baryons '92", p. 29 (1992). 

Review of Particle Physics, Phys. Rev. D 54 Part I (1996). 

O. Hanstein, D. Drechsel and L. Tiator, Phys. Lett. B 385 (1996) 45. 

R. Beck, Proc. Int. Conf. "Baryons '95", Santa Fe, (1995). 

G. Blanpied et al., preprint BNL-64382, Brookhaven 1997. 

G. Blanpied et al., Phys. Rev. Lett. 69 (1992) 1880, 

A. Sandorfi, private communication (1996), release L7a8.0. 

R. Beck et al., Phys. Rev. Lett. 78 (1997) 606. 

M. G. Olsson, Nucl. Phys. B 78 (1974) 55. 



17 



[51] 
[52] 



[53] 
[54] 
[55] 
[56] 
[57] 
[58] 

[59] 
[60] 
[61] 
[62] 
[63] 
[64] 
[65] 
[66] 



R. M. Davidson and N. C. Mukhopadhyay, Phys. Rev. D 42 (1990) 20. 
G. Hohler and A. Schulte, -kN Newsletter 7 (1992) 94. 

For a detailed discussion and interpretation of the speed of a resonant amplitude see: 

Foundations and Applications, Springer (1993). 

T. Sato and T.-S. H. Lee, Phys. Rev. C 54 (1996) 2660. 

MAMI proposal A2/2-93, Experimental check of the GDH sum rule, Mainz (1993). 
ELSA proposal to measure the GDH sum rule, Bonn (1993). 
A. Donnachie, Rep. Prog. Phys. 36 (1973) 695. 
R. Machleidt et al., Phys. Rep. 149 (1987) 1. 
P. Mergell, diploma thesis, Mainz (1995), 

P. Mergell, U.-G. Meifiner and D. Drechsel, Nucl. Phys. A 596 (1996) 367. 
M. I. Adamovich, Proc. P. N. Lebedev Phys. Inst. 71 (1976) 119. 

G. Fischer et al., Z. Phys. 253 (1972) 38. 

H. Genzel et al., Z. Phys. 268 (1974) 43. 

A. M. Bernstein et al., Phys. Rev. C 55 (1997) 1509. 

V. Bernard, N. Kaiser and U.-G. Meifiner, Z. Phys. C 70 (1996) 483. 

A2 Collaboration at MAMI, private communication (1995). 

T. Fuji et al., Phys. Rev. Lett. 28 (1972) 1672. 

G. Chiefari et al., Lett. Nuovo Cim. 13 (1975) 129. 



A. Bohm, Quantum Mechanics: 



TABLE I. Compilation of the data included in our fit. The total number of data points is 1751. 



reaction 



observable 



leV] 



rmmbcr of data points 



reference 



7P ■ 



7P- 



7n — > TT p 



da l(Kl 
da jdQ. 
da /dO. 

E 

E 

T 
da/dQ 
da/dO, 
da/dfl 

E 

T 
da/dQ 
da/dfl 
da/dfl 
da /dO, 
da/dfl 



162-220 
270-420 
247-452 
270-420 
211-267 
220-425 
160-180 
270-420 
210-425 
270-420 
272-398 
190-420 
190-230 
230-420 
250-420 
194-270 



152 

176 

61 

176 

36 

107 

224 

107 

82 

107 

28 

189 

20 

78 

158 

50 



lie] 
|i^] 
lie] 
|i^] 

|2C] 
|1E] 

lie] 

1121 

lie] 

|2C| 

\B 

121] 

\2c] 

12^1 
\2l] 



TABLE II. Vector meson coupling constants from our fit compared to the values used for the Bonn potential ]|57[ and a 
dispersion theoretical analysis of nucleon form factors |5q|. 







9^ 


9j 


9^ 


gZ 


this work 


4.85 


15.69 


6.78 


-1.67 


Ref. [g 




3.24 


19.81 


15.85 





Ref. 




1.99 


12.42 


20.86 


-3.41 
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TABLE III. The S wave amplitude £"0+ at threshold in units of 10 ^ /m^. 





7P Tv^n 


'yn ^ n p 


7P — + TT^p 




7n ^ TT n 




27.5 


-32.0 


-2.4 


0.4 


ChPT 


28.2 ± .6 


-32.7 ± .6 


-1.16 


2.13 


this work 


28.4±.2 


-31.9±.2 


-1.22±.16 


1.19±.16 


experiment 


28.3 ± .2 H 


-31.8 ± .2 Hi 


-1.31±.08 (liHI 





TABLE IV. Result of the speed plot 



Mr (MeV) Fr (MeV) 




r(10-V"^-) * (deg) 


Mi+ 1212 ±1 99 ± 2 
Ei+ 1211 ±1 102 ±2 




21.16 -27.5 
1.23 -154.7 


TABLE V. The ratio Rem and the helicity amplitudes Ai and ^3 

2 2 

for the A(1232) a total width of 113 MeV has been assumed. 


for the 7A'" — > 


A transition. In determining Ai and A3 , 

2 2 


Rem (%) Ai 




As {IQ-yVGeV) 


global fit -2.54±0.10 
local fit -2.33±0.17 


-131.2 
-129.4±1.3 


-252.2 
-246.6±1.3 



19 




20 




21 





22 




23 




24 




25 




26 




27 




28 




29 




30 




31 



-"fc 

E 



ro 
I 

O 



s -1 



-2 



—I 1 1 r- 



—I 1 1 r- 



—I 1 1 r- 



—I 1 1 r- 




100 



200 



300 
(MeV) 



400 



500 




100 



200 



300 
E^ (MeV) 



400 



500 



10 



to 

I 

O 



^ -10 



o 
LxJ 



-30 



—I 1 1 p- 



— I 1 1 1— 



—I 1 1 p- 



— I 1 1 1— 




_l I I I I I I I I I I I I I I I I I l_ 



100 



200 



400 



500 



300 
E^ (MeV) 

FIG. 13. Real (boxes) and imaginary (diamonds) parts of the isospin components of Eo+ obtained with our local fit. The 
solid lines are the corresponding results of the global fit, the dashed lines indicate the error band. 
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FIG. 16. Isospin components of Mi+. Symbols as in Fig. |l^. 
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FIG. 17. Decomposition of the isospini channel of the amplitude Eq+ into components for reactions off the proton and the 
neutron respectively. Symbols as in Fig. [13. 
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FIG. 19. Decomposition of the isospini channel of the amplitude £1+ into components for reactions off the proton and the 
neutron respectively. Symbols as in Fig. [13. 
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FIG. 20. Decomposition of the isospin i channel of the amphtude Mi+ into components for reactions off the proton and the 
neutron respectively. Symbols as in Fig. [13. 
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FIG. 21. DifFerential cross sections for n production at low energies. The results of our global fit compared to the TRIUMF 
data for jn — ^ 7r~p. The data points for Eo+ are also from Ref. [ po[ . 
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FIG. 23. The amplitude Eo+{pTr'') in comparison with the recent analyses of data from Mainz p^ , |6^ , o p en circles, and 
Saskatoon diamonds. The open squares are the result of an analysis of an earlier experiment at Mainz |^,^. The dashed 
line shows our prediction for Im_Eo+ , the solid line for Re_Bo+ respectively. The unitarity parameter j3 discussed in describes 
the rise of ImiJo+ at tt'^ threshold via Im_Eo+ ~ iP<l+- In our calculation (3 — 4.3, to be compared with the values 3.76 ± 0.11 
corresponding to the unitary fit of Ref. and the ChPT value 2.78 all in units of 10~^/m^+. 
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FIG. 24. The ratio E^' /MH' as a function of Ej. Solid line: global fit, diamonds: local fit, circles: analysis of Krahn 
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FIG. 26. Left: The speed of the multipoles Mj^ and E-^^ (solid line) compared to the case of an ideal resonance pole 
(dotted line). Right: The speed vectors for these multipoles in the complex plane (crosses) compared to the case of an ideal 
resonance pole (squares). 
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FIG. 28. Prediction of our calculation for the total cross sections (solid lines). The contributions of the leading S and P 
wave multipoles are also shown: Mi+ (dashed lines), _Eo+ (dashed-dotted lines), Mi_ (long dashed lines), Ei+ (dotted lines). 
The data are from Mainz for vr"*" and 7r°, and from Tokyo and Frascati j66| for n~ . 
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FIG. 29. Our prediction for the difference cri — as on the proton (solid line) and the neutron (dashed-dotted line) and the 
corresponding contributions to the GDH integral, which are evaluated as functions of the upper limit of integration. We also 
indicate the predictions of the sum rule. 
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